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Abstract 

We study the quantization of the noncommutative selfdual 3 model in 4 di- 
mensions, by mapping it to a Kontsevich model. The model is shown to 
be renormalizable, provided one additional counterterm is included compared 
to the 2-dimensional case [1] which can be interpreted as divergent shift of 
the field 6. The known results for the Kontsevich model allow to obtain the 
<3 genus expansion of the free energy and of any n-point function, which is finite 

for each genus after renormalization. No coupling constant or wavefunction 
renormalization is required. A critical coupling is determined, beyond which 
the model is unstable. This provides a nontrivial interacting NC field theory 
^ , in 4 dimensions. 
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1 Introduction 

This paper is an extension of our previous work [1] on the noncommutative (NC) Eu- 
clidean selfdual <p z model. In [1] we considered the 2-dimensional case, and showed 
that this model can be renormalized and essentially solved using matrix model tech- 
niques. This was achieved by mapping it to the Kontsevich model. 

The map from the selfdual NC 3 model to the Kontsevich model exists for any 
even dimension, however the eigenvalues and their degeneracy of the corresponding 
Kontsevich model depend on the dimension of the underlying space. Therefore the 
properties of the model and in particular its renormalizability must be studied sep- 
arately in different dimensions. In the present paper we elaborate the 4-dimensional 
case. Generalizing [1], we prove renormalizability and obtain closed expressions for 
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the genus expansion of the free energy F, based on results of [2,3] for the Kontsevich 
model. The general n-point functions can then be computed by taking derivatives 
of F, in an explicit way for diagonal entries <pu and somewhat implicitly for the 
general case. As an example, we work out the genus contributions for the 1- and 
2-point functions. 

It turns out that as in the 2-dimensional case, the required renormalization is 
determined by the genus sector only, and can be computed explicitly. We show that 
all contributions in a genus expansion of any n-point function are finite and well- 
defined after renormalization, for finite coupling. A linear (tadpole) counterterm 
must be introduced which is linearly divergent, while mass and a further counterterm 
are logarithmically divergent. The coupling constant does not run as expected. This 
implies but is stronger than perturbative renormalization. We thus obtain a fully 
renormalized and essentially solvable NC field theory with nontrivial interaction, by 
applying the rich structure of the Kontsevich model related to integrable models 
(KdV flows) and Virasoro constraints [2-4]. The model is free of UV/IR diseases 
due to the confining oscillator potential introduced in [5-8]. 

As in [1], these results are obtained starting with purely imaginary coupling 
constants iX, but allow analytic continuation to real coupling as long as |A| is small 
enough. If A is larger than some critical value, the model becomes unstable. We 
determine the corresponding critical coupling, which is interpreted as instability 
induced by the finite potential barrier. 

This paper is largely parallel to our previous work [1] on the 2-dimensional case. 
Therefore we will be short in certain issues which have already been discussed there, 
and which apply without change. Nevertheless, the present paper is essentially self- 
contained. In section 121 we define the 3 model under consideration, and rewrite 
it as Kontsevich model. There is one additional counterterm compared to [1]. We 
then briefly recall the most important facts about the Kontsevich model in section 
H3 Renormalization and finiteness are established in section 14.11 which is the main 
result of this paper. The 2-point function at genus is worked out explicitly in 
section I4.2| and the computation of the general n-point function is discussed in 
section 14.31 The critical point is determined in section 14.51 We conclude with a 
discussion and outlook. 



2 The noncommutative model 

We consider the 4-dimensional scalar noncommutative 3 model, with an additional 
oscillator-type potential in order to avoid the problem of UV/IR mixing. The model 
is defined by the action 

r 1 2 ■ \ 

S= ^0 2 + ft 2 (^)(^) + ?-</> 3 (1) 
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on the 4- dimensional quantum plane, which is generated by self-adjoint operators 1 
Xi satisfying the canonical commutation relations 

(2) 

i,j = 1,2,3,4. We also introduce 

Xi O^j Xjj Xj] ^ ji (^) 

assuming that 0y is nondegenerate. The dynamical object is the scalar field <fi = <fy, 
which is a self-adjoint operator acting on the representation space TC of the algebra 
P)). The term Q 2 (xi<j))(xi<j)) is included following [5-8], making the model covariant 
under Langmann-Szabo duality, and taking care of the UV/IR mixing. We choose to 
write the action with an imaginary coupling iX, assuming A to be real. The reason 
is that for real coupling A' = iX, the potential would be unbounded from above 
and below, and the quantization would seem ill-defined. We will see however that 
the quantization is completely well-defined for imaginary iX, and allows analytic 
continuation to real A' = iX in a certain sense which will be made precise below. 
Therefore we accept for now that the action S is not necessarily real. 

Using the commutation relations (J2j), the derivatives di can be written as inner 
derivatives dif = —i[xi, /]. Therefore the action can be written as 

/2 ^ 
-(xi(f)Xi(j) - XiXi4>4>) + Vt 2 xicpXicj) + —<p 2 + — 3 (4) 

using the cyclic property of the integral. For the "self-dual" point Q = 1, this action 
simplifies further to 

S = J {m + ^)0 2 + | 3 = Tr(ij0 2 + ! 3 ). (5) 

Here we replaced the integral by J = (27i6) 2 Tr, and introduce 

J = 2(27r^ 2 QTx J x i + ^), A = (2^) 2 A. (6) 

i 

We assume that 9ij has the canonical form 9±2 = —#21 ='■ 9 — #34 = — #43. Then 
J is essentially the Hamiltonian of a 2-dimensional quantum mechanical harmonic 
oscillator, which in the usual basis of eigenstates diagonalizes as 

u 2 9 

Jim, n 2 ) = 8n 2 9 (m + n 2 + 1 + ^-)K n 2 ), rn G {0, 1, 2, ...}. (7) 

: We ignore possible operator-technical subtleties here, since the model will be regularized using 
a cutoff TV below. 
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To simplify the notation, we will use the convention 



n = (ni,n 2 ), n = n x + n 2 (8) 

throughout this paper, keeping in mind that n denotes a double-index. 

In order to quantize the theory, we need to include a linear counterterm —Tr(i\)a q 
to the action (the explicit factor iX is inserted to keep most quantities real), and - 
as opposed to the 2-dimensional case [1] - we must also allow for a divergent shift 

-> + iXc (9) 

of the field 0. These counterterms are necessary to ensure that the local minimum 
of the cubic potential remains at the origin after quantization. The latter shift 
implies in particular that the linear counterterm picks up a contribution — Tr(iX)(a+ 
cJ)4> from the quadratic term. Therefore the linear term should be replaced by 
—Tr(iX)A(f) where 

A = a + cJ, (10) 

while the other effects of this shift 0^0 + iXc can be absorbed by a redefinition of 
the coupling constants (which we do not keep track of). We are thus led to consider 
the action 

involving the constants iX, a, c and fi 2 . The additional constant terms in (|TT|) are 
introduced for later convenience. By suitable shifts in the field 0, one can now either 
eliminate the linear term or the quadratic term in the action, 

S = Tr ( - -i- M 2 + l 4 3 ) = Tr ( -MX 2 + ^ X 3 - — i-M 3 ) (12) 
V 2iX v 3! v J V2 3! 3(iA) 2 / v ; 

where 2 

= + ^J = X + ^M (13) 
iX tX 

and 



M = v/J 2 + 2{iX) 2 A = ^.P + 2{iX) 2 a-(iX) A c 2 (14) 
J = J+{iXfc. (15) 

This has precisely the form of the Kontsevich model [3]. Coupling the field linearly 
to the source will be very useful for computing correlation functions. In particular, 
this can be compared directly with the 2-dimensional case considered in [1], however 
Jk plays now the role of J^, and a further constant c has been introduced. 

2 for the quantization, the integral for the diagonal elements is then defined via analytical con- 
tinuation, and the off-diagonal elements remain hermitian since J is diagonal. 
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2.1 Quantization and equations of motion 

The quantization of the model (jll)) resp. (|12J) is defined by an integral over all 
Hermitian iV 2 x N 2 matrices 0, where iV serves as a UV cutoff. The partition 
function is defined as 

Z{M) = J D0exp(-Tr(-^-M 2 0+^0 3 )) = e F ( M ), (16) 

which is a function of the eigenvalues of M resp. J. Since N is finite, we can freely 
switch between the various parametrizations ()11|). ()12|) involving M, J, 0, or 0. 
Correlators or "n-point functions" are defined through 



njn) = Z I D< P ^V(-S)<Pi lh -4i n jm (17) 



keeping in mind that each i n denotes a double- index (JBJ). The (degenerate) spectrum 
of J resp. M is given by (jZJ) resp. (|14)1. The nontrivial task is to show that all n- 
point functions have a well-defined and hopefully nontrivial limit N — > oo, so that 
the "low-energy physics" is well-defined and independent of the cutoff. 

Using the symmetry Z(M) = Z{JJ- l MU) for U G U(N 2 ) } we can assume that 
M is diagonalized with (ordered) eigenvalues m;. In fact, Z depends only on the 
eigenvalues of M 2 resp. J 2 . There is a residual U(l) x (7(2) x U(3) x ... symmetry, 
reflecting the degeneracy of J. This implies certain obvious "index conservation 
laws", such as (<f> kl ) = 5 M (4> n ) etc. 

In order to have a well-defined limit iV — > oo, we must require in particular that 
the 2-point function (4>ij4>ki) and also the one-point function (<frki) have a well-defined 
limit. We therefore impose the renormalization conditions 

Moo) = i^TT' (18) 

(0oo) = (19) 

together with the requirement that the index dependence at least of the 2-point 
function is nontrivial. This will uniquely determine the renormalization of a, c and 
/i 2 , while A receives only finite quantum corrections which will not be computed 
here. 



Quantum equations of motion and correlators. Noting that the field cou- 
ples linearly to M 2 resp. J 2 in (fT7)|) . one can compute insertions of a diagonal factor 
0^^ in a correlator by acting with the derivative operator 2i\-^ on Z resp. F . More 

general non-diagonal insertions <pki can also be obtained in principle, as discussed 
in section 14 .HI However, using various standard manipulations of the path integral 
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()16j) one can derive directly a number of nontrivial identities for the n-point func- 
tions. Since their derivation as given in [1] is independent of the eigenvalues of J, 
we simply write down the most important identities here, noting that J must be 
used instead of J. In particular, one finds for the propagator 

(<frki4>ik) = — 2 (<fec - 4>u) (20) 

m k ~ m l 

for k ^ / (no sum), where m k denotes the eigenvalues of M. Recalling that <fi = 
4> + J, this gives 

(<f>u<Pik) = ~ 2 ~ + - 2 2zA ~ 2 {<pkk - <t>u) (21) 
Jk + J{ J k — J\ 

noting that Jf — Jf = ml — mf using ()14j) . The first term is the free contribution, 
and the second the quantum correction. This "only" requires the 1-point functions 



((f> kk ) = 2i\—^ In Z(m) = —J k + ((j) kk ) (22) 
ami %\ 



which can be obtained from the Kontsevich model, as we will show in detail. Fur- 
thermore, one can derive 



) - (2zA) V (23 ) 



\o — Vnfcfc/ \ / y_> 2 

These manipulations can be generalized [1]. In particular, certain 3-point func- 
tions can be expressed exactly in terms of 1- and 2-point functions, etc. We will 
not repeat these considerations here. Instead, the finiteness of general correlation 
functions will be established directly, by showing that the appropriate derivatives 
of the (connected) generating function F(J) = lnZ(J) are finite and well-defined 
after renormalization. This could also be used to demonstrate that the model is not 
free, in contrast to [9]. We will not bother to elaborate this, since the model will be 
renormalized for finite coupling, and the lowest nontrivial term in an expansion in 
A is manifestly finite. 

3 Some useful facts for the Kontsevich model 

The Kontsevich model is defined by 



* Km \M) = e F '^ = ^—^ , r ; 7 >- (24) 



Kon i 



J dX exp j 






'} 


J dX exp | 


-Tr (Mf.) 


} 
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where M is a given hermitian iV 2 x A^ 2 matrix, and the integral is over Hermitian 
iV 2 x N 2 matrices X. This model has been introduced by Kontsevich [3] as a combi- 
natorial way of computing certain topological quantities (intersection numbers) on 
moduli spaces of Riemann surfaces with punctures, which in turn were related to 
the partition function of the general one-matrix model by Witten [10]. It turns out 
to have an extremely rich structure related to integrable models (KdV flows) and 
Virasoro constraints, and was studied using a variety of techniques. For our purpose, 
the most important results are those of [2-4] which provide explicit expressions for 
the genus expansion of the free energy. Note that A can be introduced via 



Z Kont {M) 



f dX exp j 






j> f dX exp j 


[Tr( 




f dX exp j 


{-Tr(Mfl)] 


\ 


JdX 


exp j 




\ 



_ (25) 

where X = — A 1//3 X, M = A 2//3 M, which allows to obtain the analytic continuation 
in A. 

The matrix integral in (|2*3J) and its large N limit can be defined rigorously in 
terms of its asymptotic series. This involves in a crucial way the variables [3] 

t r := -(2r + 1)!! 0ar+i, ®r ■= -TrM~ r . (26) 

r 

One can then rigorously define [2, 3] the large iV limit of the partition function 
Z Kont (M), which turns out to be a function of these new variables only, Z Kont {M) = 
Z Kont (8i). More precisely, each order in perturbation theory is a polynomial in 
finitely many of the variables t r , and becomes independent of N for N large enough. 
This provides a rigorous definition of the </> 3 model. For more details we refer to [1-3]. 

Without renormalization (i.e. for finite or zero a), 6 r is linearly divergent for 
r = 1, logarithmically divergent for r = 2, and finite for r > 3. This is a first 
indication that the model requires renormalization. 

A further crucial fact is the existence of a 



Genus expansion. As usual for matrix models, one can consider the genus ex- 
pansion 

i^Kont pKont ^ pKont (27) 

by drawing the Feynman diagrams on a suitable Riemann surface. In principle, this 
genus expansion can be obtained as a expansion by introducing an explicit factor 
N in the action, so that the action takes the form 

S'= TrN(~M' 2 <P' + ^<P' 3 )- (28) 
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However, it was shown in [2] that the F^ ont can also be computed using the KdV 
equations and the Virasoro constraints, which allows to find closed expressions for 
given g. It is useful to use the following set of variables: 

v 

I*(uo,*i) = yitfc4*^T ( 29 ) 

p>0 ^ 

where uq is given by the solution of the implicit equation 

uq = I (u ,ti). (30) 
We note that Ik can be resummed as 

(31) 

in particular 

—^ 7 w=wr Ia - (32) 

These variables turn out to be more useful for our purpose than the t r , since the 
quantities m 2 — 2uq will be finite in the renormalized model, while the t r are not. 
Using the KdV equations, [2] find the following explicit formulas: 

w Kont _ U \^ ^0 l k 1 \^ U Q \ ^ l _a_ l _b_ / q q\ 

~ 6 2^ k + 2 k\ + 2^k+l L> a\b\ 1 ' 

k k a+b=k 

Ff ont = — ln^— , (34) 

1 24 1-V V ; 



TpKont 
b 2 



5760 



+ 29 , /3/2 - +28- J| 



'1-/03 (1-/0* (1-/O 5 



(35) 



2fc+l 



etc. All Ff ont with # > 2 are given by finite sums of polynomials in /fc/(l — It)' 3 , 
the number of which is p(3g — 3) with p(n) being the number of partitions of n. 

An alternative form of F^ ont can be obtained by solving directly the "master- 
equation" (|2~H|) at genus 0. This leads to [4] 



Ft nt = ^E^-lE^-^) 3/2 -«oEK 2 - 2 «o) 1/2 

it i 

1 Vln ( (m 2 -2n )V 2 + (mg-2n )V 2 l 
2 f^ I rrti + m fe ' 
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which is equivalent to (JHHj) but more useful in our context. All sums over double- 
indices (|SJ) here and in the following are to be interpreted as 



N-l 



E = E 



i ii,i2=0 

truncating the harmonic oscillators in (J7J). The parameter uq is given by the im- 
plicit equation (}3*2*j) . This constraint can alternatively be obtained by considering 
F Kont {fhi] uq) with Mo as an independent variable, since its equation of motion 

9 -F Kont (m l ;u Q ) = l(u -I ) 2 = (37) 



du 2 

reproduces the constraint. After renormalization, all sums will be convergent for 
the physical observables as iV — > oo. 

4 Applying Kontsevich to the </> 3 model 

We need 

Z = Z Kont [M]Z free [M] exp(^y i TrM 3 ) (38) 



where 



Z^[M] = e F f- = J dXexp (-Tr (^f\ j = J] -±= J] 



mi f -"r mi + m 

% K3 



(39) 



up to irrelevant constants, so that 

-F/ree = -- ^ h^m, + rhj) (+const). (40) 
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Therefore 



F Kont + F free + ^~ 2 TtM 3 



i i 

+Y -^E ln (v / ^ 2 -2«o+ y/rh 2 k -2uo). (41) 

In the present case, the eigenvalues are given by (|2*5Jl . (JT3j) 

K = A- 2 / 3 v /j? + 2(iA) 2 A, (42) 



and the model will be ill-defined without renormalization since Uq is linearly diver- 
gent. However, note that only the combinations a/ rhf — 2uq enter in 1)31)) and ([41 p. 
which can be rewritten using (|T4j) as 



2m = \- 2/3 yfj[+2b (43) 
where 

6 = ( 2 A) 2 a - A 4/3 M - (iA) V/2. (44) 

The point is that both J and 6 will be finite after renormalization, rendering the 
model well-defined. 

The model and therefore the free energy F depends a priori on the parameters 
a, c, A, and Jk resp. which in turn contains /i 2 . In the form (|41|) . new parameters 
6 and iio have appeared, which are implicitly determined by (jBj) and the constraint 

(6 - (,A) 2 a + (.A)V/2) = - £ . (45) 



Eliminating m by (JHJ), the genus contribution to the partition function (|4~Tj) takes 
the form 

Fo = ln^ = ^^v(^%(a-^-(,A) 2 c 2 /2)^/^ 

i i 

i,k ' 

(46) 

We consider F = F(J 2 ) as a function of (the eigenvalues of) J 2 from now on, 
with further parameters a, c, A, while 6 is implicitly determined by (jJSJ). Since 



only enters through the combination y J 2 + 26, we note that the eigenvalues can be 
analytically continued as long as this square-root is well-defined. 

We can now compute various n-point functions by taking partial derivatives of 
F = Y2 g Fg (where F g = F^ ont for g > 1) w.r.t. J|, as indicated in section I2~T1 For 

the "diagonal" n-point functions (4>u---<f)kk) > this amounts to varying the eigenvalues 
of jf . In doing so, we must be careful to keep the parameters a, c constant since 
they determine the model, and note that b depends implicitly on J| through the 
constraint (|43j) . 

Some of these computations simplify by the following observation: The constraint 
(|45|) for b arises automatically through the e.o.m as in (|37|h noting that 

%-MJi) = l^Fo(Jh b) + |-Fo( J?; 6)4-6 = (# 6) (47) 



<9J 2 v w <9J 2 v *' ' <% uv l ' 'a/ 2 <9J 2 
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using 



I F °&> b) = -i ((^ - a + (zA) V/2) + £ )' = 0. (48) 



J, 2 + 26 



Thus for derivatives of order < 2 w.r.t. J|, we can simply ignore b and treat 
it as independent variable, since the omitted terms ()48|) vanish anyway once the 
constraint is imposed. 

4.1 Renormalization and finiteness 

The 1-point function We can now determine the required renormalization of a 
and c, by considering the one-point function. Using (jiSJ). and P5|). the genus 
zero contribution is 

(<M 9 =o = 2a-^F (J 2 ) 

•^+^ + -7^= ((a - ^ - (iA)V/2) 
zA * / J 2 + 26 V ^ 



A' 



% sfJl + 2b + ^r- + 2b' dJ l m 

t- ,- lX) - (49) 



and therefore 



^kk)g=0 — (4>kk)g=0 — 



L(y/fi k + 2b-J k ) + (iX)c + Y^ 



(iX) 



Jl + 2b\lJ] + 2b+(J] + 2b) 

(50) 



which must be finite and well-defined in the limit N — > oo. 

As opposed to the 2-dimensional case, we now have to face a logarithmic di- 
vergence in the sum on the rhs. We note that resp. depends only on the 
combination 

k:=h + k 2 (51) 
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(recall that k, j, ... etc. are 2-component indices (jHJ))- In analogy to the usual strategy 
in renormalization, we consider the Taylor-expansion of 



fit) := E r~ r= = f(0) + (k)f(0) + ... 

j ^J 2 + 2by/Jf + 2b+ (J 2 + 26) 
=: f(0) + f R (k) (52) 
in k, where (keeping only the divergent part) 

1 1 f N 1 

f(0) ~ > — h const ~ - — — — ■ / dxidxo -, ttt + const 

K J V + 26 (8vr 2 fl) 2 Jo (x! + a; 2 + const) 2 

ln(iV) (53) 



j 3 
1 



(8vr 2 #) 2 

up to finite corrections. We anticipate here that Jk (jBTj) and 6 (J5Bj) will be finite 
after renormalization. Only /(0) is divergent in (|52p. while all derivative terms /'(0) 
etc. and in particular are finite and well-defined as iV — > oo. 

This leads to a clear candidate for the appropriate scaling of the bare constants 
a, c, fi 2 and A: the log-divergence in (jBTIjl should be absorbed by c, and /x 2 should be 
determined such that Jk is finite. This allows to have finite a coupling constant A, 
which does not require renormalization. Keeping b also finite determines a through 
the constraint (|45jl. rendering the one-point function (jBTIj) finite and well-defined. 
Thus we set 

c = -/(0) + c' « -T^p In(iV) + c' (54) 

where d is a free finite parameter. Furthermore using 

J k = J k + {iXfc = 8ir 2 6(k + 1) + (4vr 2 ^V + (iA) 2 c) (55) 

we choose according to the above discussion 

(*A) 2 2 _ (z\) 2 f(0) - (zA)V 2 
^ " 47r 2 ^ C+ ^" 47r 2 ^ 2 

- ln(iV) (56) 



256vr 6 fl 4 

where /z 2 ^ > is finite, and a free parameter of the model. This leads to 

J k = 87r 2 e(k + 1 + ^-), (57) 
which is finite and independent of N. 
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The parameter b is determined by the renormalization conditions (0oo) — (0oo) — 
= (fT§|) . At genus 0, this amounts to 

J 2 + 2b - J + (*A)V = (58) 

using the above definitions. This has indeed a solution with finite b as long as 
| (tA j c | < 1, determining b as a function of d. Note that 6 = 0(A) 2 is always 
analytic in A, starting at second order. One particular solution is 

c ' = =>■ b = (59) 

where the formulas take a very simple form. 

Since J and 6 are now finite and independent of N, we obtain the renormalized 
one-point function 

(<f>kk) 9 =o = ^(\fj[+2b - J fc ) + (iA)c' + (i\)f R (k), (60) 

which is finite and well-defined as N — > oo. Note that there is one additional free 
parameter compared to the 2-dimensional case, given by d . For the simplest case 
d — b — 0, this simplifies as 

(<i>kk) g =o = (W^ = (<a)Et(ttt-ttt) (61) 



which coincides with the one-loop result (JHEj); we will comment on this fact later. 
Finally, a is determined by the constraint (|45|) . 

(^-a + (zA) 2 c 2 /2) = -J2 ]= — — =: (62) 



J 2 + 26 

1 /- v . I I 



dxidx 2 t r ~ — A^lnA^. 



8tt 2 6 J (xi + x 2 + 1) 8tt 2 6 

Therefore 

a = ^( & ) + + (^) 2 c 2 /2 = (7(0) + (*A) 2 c 2 /2 + (finite). (63) 
(zA) 2 

Thus we can trade a for the implicit parameter b, and interpret b as parametrization 
of a which is determined by the renormalization conditions (0oo) — 0. 

To summarize, the one-point function is renormalized by requiring the bare pa- 
rameters c, fi and a to scale as in (p>4"|) . (|5T)j). and (jHSJ- This leaves 3 independent 
finite parameters A, d and fiR of the model. A fourth free parameter b could be 
introduced by relaxing the condition (0oo) = 0. 
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Diagonal n-point functions By taking higher derivatives of F w.r.t. J 2 , we 
obtain the genus contribution to the connected part of the n-point functions for 
diagonal entries (<pi 1 i 1 ---<pi n i n ) ■ Note that the (infinite) shift 4> = 4>+j^ = 4>+j^ — (i\)c 
drops out from the connected n-point function for n > 2, therefore for n > 2 these 
coincide with (<fri 1 i 1 ---4>i n i„) and should thus be finite. 

Start with the genus contribution. To compute higher derivatives of F w.r.t. 
J 2 , we also must take into account the implicit dependence of 6 on J 2 . We recall 
that (}4~5|) can be written as 

| F„( J? ; b) = -\ (A + a + AV/2 - A-/ 3 J )" (64) 



which vanishes through the constraint (|45p. Considering this as a function of the 



independent variables J| and 6, we will also encounter 



^Fo(J?;b) = -(y 2 +a + X 2 c 2 /2-\- 2 ^I yi-\^P (J 2 ,b)), 
82 F (J 2 ;b) = (A + a + AV/2-A- 2 / 3 / )(A-/ 3 A /o( J 2 , 6)) 



dJ 2 db y 1 ' VA 2 ' V v 94 

which still vanishes due to the constraint, while e.g. 

dJ 2 dJ 2 db v dJ 2 r dJ 2 ' 

-(^ + a + AV/2-A-^/„)(A-^ 53 | 3f /„) 

is non-vanishing, and similarly for higher derivatives. We will show below that all 
these terms are finite as N — > 00 provided A is small enough. 
Consider the quantities involved in more detail. First, all the 

Wl !) = -«,- DHA^ £ (J?+ ^ + , (65) 

are finite (i.e. convergent as TV — > 00) for p > 1. In particular, we note that < 1 
provided the coupling A is small enough. Furthermore, the constraint (|43|) implies 

-t& = -_* A 2 fV^ )-Lb, (66) 

dJ 2 2(J 2 + 26)3/ 2 V If (J2 + 26)3/ 2 ; 5 J| 

hence 



«9J| 2(J| + 26)3/2 i_ 7l 
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which is finite for small |A| since we can assume < 1 as just shown. Finally, we 
need 

y p {ib) = -\-* /3 i P+1 , 

d (2p+l)!!AW)/ 3 

-=-J p (J,6) = — — — i — . (68) 



9J I 2(J| + 26r 



2 



which are all finite as iV — > oo provided A is small enough, for any p > 0. The 
same is obviously true for all higher derivatives. Furthermore, we see explicitly that 
any higher derivative of the r.h.s. in w.r.t. J| or b is also manifestly finite. 
Combining all this, we have shown that all relevant derivatives of Fq are finite and 
have a well-defined limit N — > oo, provided A is small enough, the constraint (|%H|) is 
imposed, and the parameters a, c, /z 2 are renormalized as in (jHljl . (JHHj) . and (jHHJ)- This 
establishes that the genus contribution to the correlation functions for diagonal 
entries ((j>kk--- ( t>ii) is finite and well-defined. 



Higher genus contributions. It is easy to see that the higher genus contribu- 
tions are also finite under the same conditions. Indeed, in view of the structure of 
the higher genus contributions F g stated below as found by [2], this follows 
immediately form the above considerations. Thus we have established 

Theorem 1 All derivatives of F g w.r.t. J\ forg > as well as all F g for g > 1 
are finite and have a well-defined limit N — > oo ; provided A is small enough, 
the constraint fl^<5| ) is imposed, and the parameters a, c, /i 2 are renormalized as 
m fl| ; and (Ef . 

The precise condition for A is that < 1. The limiting case I\ — 1 will be identified 
below as an instability of the model. 

Since the connected n-point functions are given by the derivatives of F — 
J2 g >o Fg w - r -t- Jki ^ ms implies that all contributions in a genus expansion of the 
correlation functions for diagonal entries (4>kk---4'ii) are finite and well-defined. The 
general non-diagonal correlation functions are discussed in section [Q| and also turn 
out to be finite for arbitrary genus g under the same conditions. Putting these re- 
sults together we have established renormalizability of the model to all orders in a 
genus expansion, i.e. 

Theorem 2 The (connected) genus g contribution to any given n-point func- 
tion is finite and has a well-defined limit N — > oo for all g, under the above 
conditions. 

Moreover, they can in principle be computed explicitly using the above formulas. 
In particular, since any contribution to F g has order at least A 4s_2 , this implies 
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renormalizability of the perturbative expansion to any order in A. This is certainly 
expected in view of the results in [5-7]. 

Note in particular that we did not have to specify whether A is real, or iX is 
real, etc. Rather, all genus g contributions are analytic in A provided |A| is small 
enough such that < 1 holds. Therefore we have obtained a definition of the NC 
<j) 3 model also for real coupling A, under the above conditions. 

It is also worth pointing out that only the genus contribution requires renor- 
malization, while all higher genus contributions are finite. This is very interesting 
because the genus contribution can be obtained by various techniques in more 
general models, see also [11,12]. This is due to the presence of the oscillator-like 
potential in the action, which suppress the UV/IR mixing originating from higher 
genus diagrams. 

The parameters fi R , A and d are the free moduli of the model, which can be 
interpreted as mass, coupling, and a further parameter which was introduced by 
renormalization. 

4.2 The 2-point function at genus 
4.2.1 (4> k i4> lk ). 

We can use (l2~Tf to obtain the genus contribution to the 2-point function ((f)ki4>ik) 
for k ^ I. Using (|60|) we obtain 



Jjl + 2b - J J* + 2b + (zA) 2 (f R (k) - f R (V>) 
(Mik)g^ = 2-¥ * n n . (69) 



P - P 

Note that A enters also implicitly through b ([58)1 . For the simplest case 6 = 0, this 
simplifies to 

{(pkl<Plk)g=Q = ~ - + 2(tA) ~- 2 ~- 2 

Jk + J I J k — J i 

= ^-^-{i\y^-^Y.—~ — kr~ — — ( 7 °) 

Jk + Jl Jk + Jl j Jj(Jk + Jj){Jl + Jj) 

The first term is the free propagator, while the second term has somewhat stronger 
decay properties. This can be rewritten in a more suggestive way using the identity 

(4>kk-4>u} , (4>kk + 4>u) _ n Jk(4>kk) - Ji(4>u} (7U 



J 2 k - Ji ^ + Ji) 2 (J k + Ji) 2 (Jk - Ji, 
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which for b = and using (JtjTj) . (jo^j) leads to 

(6 6 )- 2 2(-:A) ^ fcfc + ^ l 4( ' A)2 1 

' "° J* + Ji (J* + (■/* + ^ Y ( Jfc + ^ + J i) 

(72) 

Comparing with the perturbative computation (|94|) . this means that the one-loop 
contribution to the propagator gives the exact result for the genus sector for the 
case 6 = 0, noting that A 2 c = 5 J using (jHfij) and (|87jl. This remarkable property can 
be traced back to (|4fi|). which is also exact at order A 2 for 6 = 0. For higher genus, 
this is no longer the case. 



4.2.2 (6 H 6 kk ). 



As a further example, consider the 2-point function (6u6 kk ) for k ^ I, which vanishes 
in the free case. To compute it from the effective action, we need in principle 



{6 a 6 kk )-(6 kk )(6 u ) = 2i\- 1 -2i\-^-(F + F 1 + ...). (73) 

dJf dJf: 

Even though this corresponds to a nonplanar diagram with external legs, it is ob- 
tained by taking derivatives of a closed genus ring diagram. Therefore we expect 
that only F will contribute, and indeed the derivatives of F\ contribute only to 
order A 4 . We need 



d d 
2i\^^2i\— ^F 
dJ? dJl 



2{i\f p _d_ 
Jl + 2b 1 d ~ J i 



( 



E 



i 



J t 2 + 26- 



J 2 



2b 



1 



= (*A) 2 r _ 

sjJl + 2byjj? + 2b 
using again (|4*K|) . Therefore to lowest order we obtain 

1 




(74) 



Jk Jl \Jk + Jl, 

in complete agreement with the perturbative computation (jHSJ- 



(75) 



4.3 General n-point functions 

Finally we show that all contributions in the genus expansion (and therefore per- 
turbative expansion) of the expectation values of any n-point functions of the form 

(6 ilh ....6 injn ) (76) 
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have a well-defined and finite limit as iV — > oo provided b is finite, which means that 
the model is fully renormalized. The argument is a generalization of the one given 
in [1], taking into account the degeneracy of J. 

In view of f{lfij) . the insertion of a factor fa can be obtained by acting with the 
derivative operator 2iA^L~ on Z(J 2 ) , resp. F g (J 2 ) for fixed genus g. Now any given 

correlation function of type (fT6|) involves only a finite set of indices i,j, ... . Thus 
taking derivatives w.r.t. Jfj amounts to considering matrices J of the form 



J 



diag(Ji, ... J fc ) + 5J kxk 








diag(J fc+ i, ...Jjv) 



(77) 



where k is chosen large enough such that all required variations are accommodated 
by the general hermitian k x k matrix 



Jkxk ■= fdiag(Ji, ...J k ) + 5Jfcxfc) (7£ 



in ([77]) . while the higher eigenvalues J k +i, ... Jjv are fixed and given by (JH7|) . Therefore 
we can restrict ourselves to this k x k matrix, which is independent of N. As was 
shown in section 14. 1[ all F g are in the limit N — > oo smooth (in fact analytic) 
symmetric functions of the first eigenvalues squared, hence of the eigenvalues of 
(Jkxk) 2 - Such a function can always be written as a smooth (analytic) function of 
some basis of symmetric polynomials in the J 2 , in particular 

F g (J 2 ,...J 2 k ) = f g (Tr(J 2 kxk ),...,Tr(J 2 kxk )). (79) 

This can be seen by approximating the analytic function F g (zi, ..,z k ) at the point 
Z{ = Jf by a totally symmetric polynomial in the Zi, which correctly reproduces the 
partial derivatives up to some order n. According to a well-known theorem, that 
polynomial can be rewritten as polynomial in the elementary symmetric polynomi- 
als, or equivalently as a polynomial in the variables s n := ^ zf, n = 1, 2, k. This 
amounts to the rhs of (|79*]). 

In the form (|79jh it is obvious that all partial derivatives of F g exist to any 

given order, and could be worked out in principle. This completes the proof that 
each genus g contribution to the general (connected) correlators (faj l ----4>i n j n ) is 
finite and convergent as iV — > oo. This implies in particular (but is stronger than) 
renormalizability of the perturbative expansion to any order in A. 



4.4 Approximation formulas for finite coupling 

In this section we derive some closed formulas which are appropriate for finite cou- 
pling A, in the large N limit. They will be needed in particular to derive the critical 
line in section l4~Hl For simplicity we only consider the case b = d = (|59jl . 
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We approximate the various sums by integrals. This gives 

^2/3 r%N fVN l 



J ° = " A / Et^I / dx dy 



A 2/3 



where 



Similarly, 



NlnN, (80) 



x n = (n H zLifi^ dx = dn . 



, 2 v- 1 A 2 p 1 



(J2)| (8^)3 7^ "(s + y+l)* 

a- , i a 2 i . . 

JTT^r « - O , _ 2m3 .on I o - ( 82 ) 



2(8**0)3 4 (a . + ( i±|M } + 1)2 2(87r 2 ^)3 tfrf + 2' 



Similarly, all J p can be approximated by elementary, convergent integrals. 
4.5 Critical line and instability. 

We have seen that for small enough coupling |A|, the free energy F = F Q + F\ + ... 
is regular and finite for any given genus in the renormalized model, since all Ik with 
k > 1 are finite provided 

However, as is manifest in the explicit formulas for F g at higher genus (J3"3j) ff-> 
there is a singularity at Ji = 1. Using ()82j) . this critical point is given for d = b = 
by 

"«' , + 2 =2» < 83 > 

This 3 is similar to the 2-dimensional case [1], and indicates that for the 3 model 
with real coupling constant A' = iX stronger than this critical coupling, the model 
becomes unstable. This is very reasonable, since the potential is unbounded, and the 
potential barrier around the local minimum becomes weaker for stronger coupling. 
Therefore this critical line could be interpreted as the point where the quantum 
fluctuations of (f> become large enough to see the global instability, so that the field 
"spills over" the potential barrier. Similar transitions for a cubic potential are known 
e.g. for the ordinary matrix models, but may also be relevant in the context of string 
field theory and tachyon condensation [13,14]. In particular, it is interesting to note 

3 Recall that this is obtained imposing the renormalization conditions (</>oo) — at genus 0. 
Note also that the rhs is indeed dimensionless. 
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that this singularity occurs simultaneously for each given genus, which suggest that 
some double-scaling limit near this critical point can be taken, again in analogy with 
the usual matrix models (for a review, see e.g. [15]). Again, such a scaling limit for 
the Kontsevich model is discussed in [2]. We leave this issue for future work. 

5 Perturbative computations 

We write the action (fTTj) as 

S = Tr(i(J0 2 + 2 J) + ^0 3 -(a)A0) 

= Tr{\<P) (G K )g 4t + ^ 3 - (zA)A0 + 1(5J0 2 + <f> 2 6J)) (84) 

where the finite (renormalized) kinetic term {Gr){'^. = ^^{J^ + Jf) defines the 
propagator 

\i\k _ (jAxk\ _ Zisk 2 _ ricfc 1/(47T 9) 

" } ~ Jf + Jf ~ i + (&0 + 2) ' (85) 

and we use again the notation ([51)1 n = n\ + n^- This corresponds to the finite 
(renormalized) matrix 

J R \n u n 2 ) = 8n 2 9(n+l + ^-)\n 1 ,n 2 ). (86) 

The remaining 

SJ\ ni ,n 2 ) = (8n 2 9 6 -^-)\ ni ,n 2 ) (87) 
is part of the counter-term, where Sfj, 2 = (/i 2 — fi 2 R ). 

1-point function The one-loop contribution to the 1-point function gives 
i , \ ^ a iA 1 v-^ 2 _,, 2 v 

- -%{-±+m^. i+k+ \ + a) +0(A2) ' (88) 

To proceed, we expand 

fta : = £ ititW = " (0) + tt) "' (0) + k «® (89) 
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where 



*'«» = -g^S tt+^+^ —aPg^ < 90 > 

and /ir(£) is a finite nontrivial function of i. 

We note in particular that the z- dependent term (£) h'(0) in f)88j) forces us to 
introduce a corresponding counterterm to the action, which we chose to be A = a+cJ 
(flUjl. As discussed in sectional this is equivalent to an infinite shift (Q of <p. Taking 
this into account we have 

(0»> = ~(-(a + cJO + (/i(0) + (i)V(0)+Mi))) +0(A 2 ), (91) 
and the condition (0oo) = implies 

MO), 

^"'<°>- (92) 

This is in complete agreement with (|54j). and also with (J53j) taking into account 
()96|) . These renormalization conditions guarantee that the one-point function (tpa) 
has a well-defined and nontrivial limit N — > oo. 

2-point function Next we compute the leading contribution to the 2-point func- 
tion (4>u<f)kk) for / 7^ k, which vanishes at tree level. The leading contribution comes 
from the nonplanar graph in figure ^ which gives 




Figure 1: one- loop contribution to (<pu<pkk) 
(</>«« = (« (<Pu) + (j^j) 2 ( 93 ) 



a + 8n 2 6c 
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(for I ^ k) indicating the symmetry factors, where the disconnected contributions 
are given by (|88j). This is in complete agreement with the result ()75|) obtained from 
the Kontsevich model approach. Note that the counterterm S J does not enter here. 

Similarly, the leading contribution to the 2-point function (4>ki4>ik) for I 7^ k, has 
the contribution indicated in figure El but now the counterterm 5 J does enter also. 




Figure 2: one-loop contribution to 



This gives the result 



?ki<Pik 



J? + Jl R 
4 



2(iA) 



■>kk 



! + 4>ll) 



E 



z'A) 2 



1 



SJi + 5J fc 



+ 0(A 4 ) (94) 



The first term is the free propagator, the next term the tadpole contributions, and 
the last them the one- loop contribution in figure 121 with counterterm 8 J. 

We have to adjust the parameters such that the result is well-defined and non- 
trivial. The last term is logarithmically divergent, 



E 



E 



(8^)^(j+^ + 2)2 



1 



ft'(O) 



logiV. 



Therefore the divergent terms 



— 87T 



SJi+8J k 
2 

5^ r 



8n 2 6 v ' (8ir 2 9) 2 
(^A) 2 g^ ^'(0) must cancel, i.e. 
1 



(95) 



(zA) 



87r 2 e 



256 tt 6 # 4 



logiV 



(96) 



in complete agreement with ()56|1. Hence the mass is log-divergent as expected, and 
no wavefunction-renormalization Z is required. 
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We note in particular that these one-loop computations already give the exact 
results for the renormalization, as found using the Kontsevich model. This reflects 
the "super-renormalizability" of the model, which is thus established rigorously. 

6 Discussion and conclusion 

We have shown that the selfdual NC 3 model in 4 dimensions can be renormalized 
and essentially solved in terms of a genus expansion, by using the Kontsevich model. 
This provides closed expressions for the free energy and certain n-point functions 
for each genus, which are finite after renormalization and valid for finite nonzero 
coupling. Remarkably, the genus contribution turns out to be exact at one loop 
in a special point b = d = of moduli space. An instability is found if the coupling 
constant reaches a critical coupling, as expected for the 3 model. 

It is very remarkable that a nontrivial 4-dimensional NC 3 field theory allows 
such a detailed analytical description. There is no commutative analog to our knowl- 
edge, which shows that the noncommutative world in some cases is more accessible 
to analytical methods than the commutative case. Furthermore, these NC 3 field 
theories in different (even) dimensions can all be mapped to the same Kontsevich 
model Z(M), for different M with different eigenvalues and degeneracies. While the 
techniques used in this paper are more-or-less restricted to the 3 interaction, it is 
worth pointing out that the renormalization is determined by the genus contri- 
bution only, which is accessible in a wider class of models; see also [11, 12] in this 
context. 

Perhaps the main gap in our treatment is the lack of control over the sum over 
all genera g. While the contributions for each genus are manifestly analytic in the 
coupling constant A, we have not shown that the sum over g converges in a suitable 
sense. This would amount to a full construction of the model. However, it seems 
very plausible that this is the case, and the sum defines an analytic function in A 
near the origin. It should be possible to establish this using the relation with the 
KdV hierarchy or the relation with topological gravity, which is beyond the scope 
of this paper. Another interesting extension would be the case Q ^ 1, which could 
be considered as a perturbation around Q = 1, rather than around Q = 0. We hope 
to study this problem in the near future. 
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